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Abstract 

This work is a comprehensive investigation of the Yang-Mills con- 
densate (YMC) dark energy (DE) model, which is extended to 
CO ■ 

include the 3-loop quantum corrections. We study its cosmic 

O 

evolution and the possibility of crossing phantom divide w = — 1, 
examine in details the Hubble parameter H, the deceleration pa- 
rameter q, the statefinder diagnosis (r, s), and the w—w' diagnosis 

00 ' 

of the model without and with interaction, and compare our re- 
sults with other DE models. Besides, by using the observational 
data of type la supernovae (SNIa), the shift parameter from cos- 
mic microwave background (CMB), and the baryon acoustic os- 
cillation (BAO) peak from large scale structures (LSS), we give 
the cosmological constraints on 3-loop YMC model. It is found 
that the model can naturally solve the coincidence problem, and 
its prediction of the afore-mentioned parameter is much closer to 
the ACDM model than other dynamics DE models; the intro- 
duction of the matter-DE interaction will make the YMC model 
deviating from the ACDM model, and will give an equation of 
state (EOF) crossing —1. Moreover, it is also found that, to fit 
the latest SNIa data alone, the ACDM model is slightly better 
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than the 3-loop YMC model; but in fitting of the combination 
of SNIa, CMB and LSS data, the 3-loop YMC model performs 
better than the ACDM model. 

PACS numbers: 95.36. +x, 98.80. Qc, 98.80.-k 
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1. Introduction 



Observations on the type la Supernova [1], CMB anisotropics [2] and large scale structure 
[3] all indicate the existence of mysterious dark energy that driving the current accelerating 
expansion of Universe. To interpret the physics of dark energy, there have been a large amount 
of models proposed. The simplest one is the ACDM model, which can fit the observations 
so far, but is plagued with the fine-tuning problem and the coincidence problem [4]. While 
the former problem exists for almost all the DE models, the latter one can be solved in the 
frame of dynamical DE models. Among them are quintessence [5], phantom [6], k-essence 
[7], quintom[8], tachyonic [9], holographic [10], agegraphic [11], ect. There are also other 
interesting models either based on the effective gravity [12], or on the Born-lnfeld quantum 
condensate [13]. In our previous works [14] a dynamic model is proposed, in which the 
renormalization-group improved effective YMC serves as the dark energy; and in a latest work 
[15], this model was extend to the 2-loop quantum corrections. Unlike the scalar models, 
our model is based on a vector-type quantum effective Yang-Mills (YM) fields and does not 
suffer from the difficulties of scalar models mentioned in Ref.[16]. The effective YMC in our 
model is a coherent boson field system at low temperatures, and the DE is viewed as the 
ground state energy of this YMC field. As is known, a quantum system of bosons with a 
conservative "charge" (such as the particle number, the electric charge, the color in QCD, 
etc) will experience Bose-Einstein condensation when the temperature is low enough or the 
charge density is high enough. This applies to many systems, such as Q-balls [17, 18], charged 
relativistic scalar bosons [19, 20], and the gluons condensate in the effective QCD models 
[21, 22]. And this is the physical origin of our Yang-Mills condensate model. From the 
viewpoint of quantum effective field theory at low temperatures, it would be desired if one can 
include high order quantum corrections as much as possible. In this work we will extend the 
YMC DE model to the 3-loop quantum corrections, focusing on its cosmic evolution and the 
issue of crossing w — — 1. 

As an important next step, one needs to confront DE models with observational data. 
In a recent work, by using the differential ages of passively evolving galaxies, Simon et al. 
gave 9 observational H(z) data points in the range < z < 1.8 [23], which have been 
used to constrain various DE models [24, 25]. Besides, the deceleration parameter q(z), 
containing the second order derivative of the scale factor a{t), is also important in confronting 
DE models [26, 27, 28]. Moreover, a new geometrical diagnosis pair (r, s), called statefinder, 
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has been introduced to distinguish DE models [29]. Since the pair (r, s) contains the third order 
derivative of a(t) and their values is expected to be available from the future SNAP observation 
[30], the statefinder diagnosis has attracted a lot of attention [31, 32, 33, 34, 35, 36, 37, 38, 
39, 40, 41]. Finally, another dynamical diagnosis w — w', consisting of the EOS and its time 
derivative, are also extensively used in the literatures [42, 43, 44, 45, 46, 47, 48, 49, 50, 51]. In 
viewing these different diagnosis methods, we shall present a comprehensive analysis for H(z), 
q(z), (r,s) and w — w', respectively, in the 3-loop YMC model. Comparing with the issues 
listed above, it is more common to test various DE models by using the observational data of 
SNIa. Especially, in a recent work, Riess et al released the up-to-data 182 gold sample of SNIa 
[52], which has been used in Refs. [53, 54, 55]. In addition, as useful complements to that 
of SNIa, the shift parameter R from CMB observations [56, 57] and the BAO peak parameter 
A from LSS [58] are also crucial to constrain various DE models. So in this work, we shall 
utilize the combination of SNIa, CMB, and LSS data, to give the cosmological constraints on 
our model. 

The organization of this paper is as follows. In section 2, we extend our previously proposed 
YMC DE model to the 3-loop quantum corrections, present its prediction of cosmic evolution, 
and explore the possibility of crossing the phantom divide w — —1 . In section 3, we study the 
Hubble parameter H, the deceleration parameter q, the statefinder diagnosis (r, s), and the 
w — w' diagnosis in 3-loop YMC model, and compare our results with other DE models. In 
Section 4, the observational data of SNIa, CMB, and LSS are employed to give the cosmological 
constraints on the 3-loop YMC model. Section 5 is a short summary. In this paper the unit 
with c = k — 1 is used. 

2. 3-loop YMC Model 

In the renormalization-group improved effective YM field theory [59, 60], the running 
coupling constant up to the 3-loop [61, 62], should have the following form 

9\F) = i 

where r = \n\F/eK 2 \, F = — \F a iliV F ailv = E 2 — B 2 plays the role of the order parameter of 
the YMC, and the parameter k is the renormalization scale with dimension of squared mass. For 
the gauge group SU(N) without fermions, b = rj = ^ ~ q.8. Here C = ^ - 1, 
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and the numerical coefficients (3q, f3±, (5 2 are given in Ref.[62]. It should be stressed that, 
although the 1-loop and 2-loop corrections are uniquely fixed, the 3-loop correction to g 2 
is renormalization-scheme dependent, so is the coefficient C [63]. Notice that the Lorentz 
invariance is reserved in the effective YM theory, because the Lagrangian is constructed out of 
the combinations of F a ilv F ailv [21, 22, 64]. For simplicity, we only discuss the case of pure 
"electric" condensate with F = E 2 (the case of including magnetic component was discussed 
in Ref. [65]). The effective Lagrangian, defined as L e ff = F/2g 2 (F), is given by 



1, 



L eff = 2 6/t e 



(y — 1) + i] In \y — 1 + S\ — rf 



An 2 \y - 1 + 5\ - \n\y - I + 5\ 



y-l + S 



(2) 



where the variable y = r + 1 = \n\F/n 2 \, and 5, a dimensionless parameter, represents the 
higher order corrections, including the type of terms, such as the term of C/r 3 in Eq.(l). In 
the bracket of Eq.(2) the 77 term is the 2-loop contribution, and the rj 2 term is the 3-loop 
one, which adds no new parameter other than the 2-loop model. Notice that the YM field 
introduced here is not the gluon fields in QCD, nor the gauge boson fields in the electro-weak 
unification. From the effective Lagrangian in Eq.(2), follow the energy density and the pressure 
of YMC DE 

! ' \{y + 1) + fifr + 2Y 2 ) - V 2 (Y 3 - 2F 4 ) 



Py = 



Py 
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bK 2 e y 



(y - 3) + 7j(yi - 2Y 2 ) - r] 2 (Y 3 + 2F 4 ) 



where 



Yi = 1^-1 + 51, Y 2 = 



Y 3 = (Y 1 - \)Y X Y 2 , Y 4 = (Y 1 - 3)Y 1 Y 2 2 . 
As a consistency check, the well-known conformal trace anomaly [21, 66] 

1 



F^ fi — Py 3p y 



2F— 

dr 



9 2 (F) 



(3) 
(4) 

(5) 
(6) 

(7) 



is satisfied up to the 3-loop. It is known that the trace anomaly occurs as a quantum effect of 
the vacuum and only violates the traceless condition of the stress tensor F^ u without violating 
the Lorentz invariance. Also the form of the stress tensor F^ of YM fields is consistent with 
homogeneity and isotropy of the Universe. The EOS for the YMC is given by 

Py 



w = 



Py 



(8) 
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If one ignores the terms of rf from Eq.(l) through Eq.(4), the 2-loop model [15] is obtained, 
and if one further sets 77 = 0, the 1-loop model [14] is recovered. 

In our model the Universe is filled with three kinds of major energy components: the dark 
energy represented by the YMC, the matter (baryons and dark matter), and the radiation 
(consisting of CMB and other massless particles). The overall cosmic expansion is determined 
by the Friedmann equation 

.a. o 87rG . . . . 

(-) = ^r(py + Pm + Pr), (9) 

where p m is the energy density of the matter, and p r is of the radiation. The dynamical 
evolutions of the three components are given by 

Py + 3-(p y + Py) = -Tp y , (10) 

(X 

a 

Pm + 3-p m = Tp y , (11) 
(X 

Pr + 3-{pr+Pr) = 0, (12) 
(X 

where p r is the radiation pressure, T is the decay rate of the YMC into matter, a parameter 
of the model. If T = 0, the YMC does not couple to the matter; if T > 0, the interaction 
term Yp y in Eqs.(10) and (11) represents the rate of energy transfer from the YMC to the 
matter. The sum of Eqs.(10), (11), and (12) guarantees that the total energy of the three 
components is still conserved. Replacing the old variables t, p m , p r , and p y with new variables 
N = lna(t), x = p m /\bK 2 , r = p r /\bK 2 , and y, and making use of the Friedmann equation 
at z — 0, the set of equations, (10) through (12), take the following form: 

dy _ 4 [y + i](Yi + Y 2 ) - V 2 (Y 3 - y 4 )] + [(y + 1) + y(Y 1 + 2Y 2 ) - V 2 (Y 3 - 2Y 4 )} 
dN~ (y + 2) +r] (Y 1 + 3Y 2 -2Y 1 2 )- V 2 (Y 3 - 3F 4 + 4(r x - 4)F 1 F 2 3 ) 

dx TCo 



dN H ( 



(13) 

(y + 1) + 7j(yi + 2Y 2 ) - V 2 (Y 3 - 2Y 4 )] - 3x, (14) 



ff.T 



where 

C = yfev[(y + 1) + ^ + 2Y 2 ) - V 2 (Y 3 - 2Y 4 )} +x + r, (16) 

Co = C( z — 0)> an d = H(z = 0). Once the parameters T and 5, as well as the 
initial conditions, are specified, the solution of this set of equations follows immediately. As 
our calculation will show, for the YMC to be a sensible model of dynamic DE, the order of 



magnitude of the decay rate T should be less than, or at most of order of the expansion rate, 
i.e., T < H . To be specific, we take the decay rate T = 0.31if and the parameter 5 = 4. 
The initial conditions for Eqs.(13) through (15) are chosen at a very high redshift Zi = 10 10 
during the Big Bang nucleosynthesis (BBN) era. To ensure the equality of radiation-matter 
occurring at a redshift z = 3454 [2], the initial radiation and matter are taken as 



Xi 



1.22 x 10 29 , n = 3.52 x 10 35 . (17) 



Besides, to ensure the BBN occurs as usual [67], the initial YMC fraction should be ~ 10% 
or less [15]; for concreteness we take the upper limit 

Vi < 74, i.e., ^<3x 1(T 2 . (18) 

Pri 

In Fig.l, We plot the dynamical evolution of p y , p m , and p r in 3-loop YMC model without 
and with interaction. For a whole range of initial = (1, 74), which corresponds to the initial 
energy fraction of YMC ^ ~ (2.5 x 1(T 35 , 3 x 10~ 2 ) ranging ~ 33 orders in magnitude, YMC 

Pri 

always has a desired tracking solutions, i.e., during the radiation era the YMC follows the 
radiation as p y oc p r oc a(t) -4 , then during the matter era it follows the matter approximately 
as p y oc p m oc a(t)~ 3 . Rather later around z ~ 0.5 it becomes dominant, and then it levels off 
and becomes asymptotically a constant for z < 0. The existence of a tracking solution can be 
analytically proved, also. For any DE model, the energy density and the pressure of DE can 
be written as 

Py = E k + V(y), p y = E k -V(y), (19) 

where E k denotes the kinetic energy, and V(y) denotes the potential energy. For our model, 
one can easily obtain V from Eqs.(3) and (4). Following the Ref.[68], we introduce a key 
function 

V"V 

A s J\Fy-' < 20 > 
whose properties determine whether tracking solutions exist. The prime means derivative with 
respect to y. Steinhardt et.al. [68] had demonstrated that, the tracking solutions exist if: (1) 
A is nearly constant, i.e., | A" 1 ^"^ A( y/ /V) |< 1; and (2) A > 1 for w < w B or A < 1 
for w B < w < (1/2)(1 +wb)- (Here w B denotes the EOS of background component) Our 
calculation shows that for the range y > 1 in YMC model, | 1^ ^ ' s a ' wa Y s satisfied; 

besides, w < w B and A > 1 also hold true in our model. Therefore, there is a tracking solution 
exists in our model. For a smaller initial value p yi , p y (t) still tracks p r (t), but for a shorter 



period correspondingly, and then approaches to the same constant. When the initial value yi 
is sufficient small, the YMC dark energy is effectively similar to the cosmological constant A. 
For the non-interaction case, the matter component retains p m oc a(t)~ 3 and always decays 
with time t. For the case of YMC decaying into matter, the matter density p m deviates from 
oc a(t)~ 3 around z ~ and becomes a constant at last. For the decay rate T = 0.31H , 
the dynamical equation for (x,y) at t — > oo has a fixed point (xf,yf) = (0.052,-0.887), 
regardless of initial conditions. Let we study the stability of this fixed point analytically. Base 
on Eqs.(14) and (13), one can obtain, 

m = f{x ^ m = 9{x > y) > (21) 

where f(x,y) and g(x,y) are the terms on the right hand side of Eqs. (14) and (13), respec- 
tively. By a standard procedure, expanding x = xj + e and y — yf + r\ (e and r) are small 
perturbations around the fixed point), and keeping up to the first order of small perturbations, 
Eq.(21) is reduced to 

Mi) =M {*y (22> 

where M is a 2 x 2 matrix, whose elements are M n = df % f f f) , M 12 = df % f y Vf) , M 21 = 
dg( ~ X Q f f^ , and M 22 = dg{ ~ X Q y v ^ , respectively. The general solution for the linear perturbations 
is of the form 

£ = C ie » lN + C 2 e^ N , r] = C 3 e» lN + C 4 e» 2N , (23) 

where C\, C 2 , C 3 , and C 4 are constants, p x and p 2 are the eigenvalues of matrix M. As 
long as iii and p: 2 are both negative, the fixed point is stable, and the solution is an 

attractor. By Calculating, we find the matrix 

( -3.16068 0.18234 \ 
~ ^ 0.20061 -2.55983 ) ' ^ > 

and its two eigenvalues are \i\ = —2.50411 and p: 2 = —3.21640, respectively, both negative. 
Thus the fixed point of this model is stable against perturbation, and the solution is an 
attractor. Therefore, the Big Rip would not happen in our model. For an initial YMC DE 
ranging ~ 33 orders in magnitude, and for both cases without and with interaction, the 
fractional densities n y0 = 0.73 and f2 m0 = 0.27 are always achieved at z — 0. In addition, 
from the theoretical point of view [68], for any DE model with an attractor-like behavior, as 
long as w < w B is satisfied, the coincidence problem can be naturally solved. Both these two 
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conditions are satisfied in our model. Therefore, the coincidence problem is naturally solved 
in 3-loop YMC model. 

In Fig. 2, we plot the evolution of EOS w(z) in 3-loop YMC model without and with 
interaction. At the early stage of Universe (high energies limit), w(z) approaches to that of 
radiation, i.e., w — > 1/3. With the expansion of Universe and the decreasing of energy scale, 
w smoothly decreases, and the YMC component transits from radiation to matter, and to DE. 
For the non-interaction case, w does not cross, only asymptotically approaches to —1. For the 
interaction case, w can smoothly cross w = — 1, as indicated by the preliminary observational 
data [69, 70, 71]. Adopting the initial y { = 74, the EOS of YMC will cross -1 around z ~ 0.6. 
For a smaller yi the crossing occurs earlier. For example, taking y^ = 72 in our model will 
make w crossing —1 occur around z ~ 1.5, as suggested in Ref.[25]. Moreover, treating 
T as a parameter, we find that a constant interaction T corresponds to a constant present 
EOS w , and a larger T yields a smaller w . For instance, T = 0.31H — > w = —1.05; 
r = 0.67# -> = -115; and T = 0.82H -> w = -1.21. 

In comparison with the 2-loop YMC model, the transition to the DE-dominant era occurs 
at z ~ 0.5 in the 3-loop model, later than z ~ 0.6 predicted by the 2-loop model [15]. Besides, 
3-loop model yields an EOS w = —1.06, which is closer to the observational constraints on 
w [69] than that of the 2-loop model [15]. In addition, the 3-loop model can also predict a 
larger age of the Universe than the 2-loop model [72]. Notice that the YMC is subdominant 
during the early stages, so the nucleosynthesis and the recombination occur as in the standard 
Big Bang cosmology. Besides, since the DE becomes dominant at very late era, the matter 
era is also long enough for the structure formation. It should be mentioned that the scale k 
can be fixed by requiring p y in Eq.(3) be equal to the dark energy density ~ 0.73p c , where p c 
is the critical density, yielding k 1 / 2 ~ 7.5h^ 2 x 10~ 3 eV (h is the Hubble parameter). At the 
moment we do not have an answer to the question why k is so small, so the fine-tuning problem 
is still present in our model. As has been shown [14], in the case of the YMC decaying into 
matter and radiation, both p m and p y will asymptotically approach to their respective constant 
values, i.e., the future of the universe is a steady state, quite similar to that of the Steady 
State model [73]. Therefore, in a sense, our model bridges between the Big Bang and the 
Steady State model. 



3. The diagnosis of H, q, r — s, and w — w' in the Model 
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In this section the decay rate is taken to be T = 0.31if as before. First, let us discuss 
the Hubble parameter H = a/a. The expansion of Universe is determined by the Friedmann 
equations 

H 2 = -g-p, (25) 

a AnG . . 
- = ^(p + 3p), (26) 

where the total energy density p = p y + p m + p r , and the total pressure p = p y + p r . The 
dynamical evolution of p y , p m , and p r have already been given in the previous section, so the 
Hubble parameter H can be easily obtained. In Fig. 3 we compare the observed expansion rate 
H(z) [23] with that predicted by ACDM model [23, 25] and by coupled 3-loop YMC model. 
The area surrounded by two dashed lines shows the 68% confidence interval [2]. It is seen that 
the coupled YMC model is quite close to the ACDM model in the range z < 1, and both two 
models approximately agree with the observations. The observed dip of H(z) around z ~ 1.5 
[23] is difficult for both models, but H(z) in our model is slightly lower than that in ACDM 
model and is closer to the dip. For the non-interaction case in the 3-loop YMC model, H(z) 
is more close to the ACDM model. Since its evolution trajectory is almost overlaps with that 
of ACDM model, we do not plot its curve here. 

Next, we turn to the deceleration parameter q(z), which is given by 

q = = ^(1 + 3tt y w + ft r ), (27) 

where Q y = p y /p and f2 r = p r /p- In deriving this expression, the Eqs.(8), (25), and (26) 
are used. In Fig. 4 we plot q(z) in 3-loop YMC model without and with interaction. In both 
cases, starting from a positive value during the matter era, q(z) decreases with the expansion 
of Universe, turns into negative around z ~ 1, and approaches to an asymptotic value q — — 1 
in future. The current value is q = —0.572 for the non-interacting case, which is denoted by 
a square dot; and is q = —0.656 for the interacting case, which is denoted by a round dot. In 
comparison, the ACDM model with f2 A = 0.73 has q = —0.595, denoted by a star symbol, 
so the non-interaction YMC model is closer to the ACDM model than the coupled YMC. 
Now, we study the statefinder diagnosis defined as [29] 

r "^' S ^w^wr (28) 

Taking time derivative of Eq.(26) and making use of Eqs.(10), (11), and (12), one obtains 

9 3 3r 

r = 1 + -fL y w(l +w)- -tt y w' + 2VL r + —£l y w (29) 
2 2 2x3 
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s = 



3QyW(l + W) — VtyW' + 



w 



(30) 



where 



dw dw dy 



(31) 



dN dy dN 



can be calculated by Taking variable y derivative of Eq.(8). The expressions of r and s in 
Eqs.(29) and (30) hold actually for a generic DE model. Since different cosmological models 
exhibit qualitatively different trajectories of evolution in the r — s plane, the statefinder is a 
useful tool to distinguish cosmological models [29]. In Fig. 5 we plot the evolution trajectories 
of statefinder in r — s plane for 3-loop YMC model without and with interaction, starting 
from the redshift z = 3. The arrows alone the curves indicate the direction of evolution. The 
overall profile of statefinder diagnosis predicted by these two model are quite similar. As the 
Universe expands, s increases to a maximum and r decreases to a minimum; after that, the 
trajectories turn a corner and approach to a final fixed point (r, s) = (1,0). This fixed point 
is same as that predicted by ACDM model [29], denoted by a star symbol on the plot. The 
current value of statefinder is (r, s) = (0.972,8.660 x 10~ 3 ) for the non-interacting case, and 
is (r, s) = (0.912,2.528 x 10~ 2 ) for the interacting case. Base on these calculated values, 
one sees that the non-interacting 3-loop YMC model is closer to the ACDM model; and the 
introduction of interaction between matter and DE cause a deviation from the ACDM model. 
We also find that a larger interaction T yields a larger deviation. For instance, T = 0.67H 
yields (r, s) = (0.851, 3.898 x 10 -2 ) at z = 0, which is further away from the value of ACDM. 
In comparison, the quietessence model [30] gives the current values (r, s) = (0.4,0.3), the 
Chaplygin gas model [31] gives (r, s) = (1.95,-0.3), and the agegraphic model [40] with 
parameter n = 2.0 give (r, s) ~ (—0.2,0.5). All the current values of (r, s) predicted by 
these three types of DE models are far away from the (1,0) given by ACDM model. The 
holographic model [41] without interaction gives (r, s) ~ (0.94,0.01), when an interaction is 
included, it gives (r, s) ~ (0.75,0.09), deviating away considerably from the ACDM model 
again. Therefore, our model is much closer to ACDM than other dynamics DE models. 

Finally, we investigate the w — w' diagnosis, defined in the Eqs.(8) and (31). In Fig. 6 we 
plot the evolution trajectories of w — w' for 3-loop YMC model without and with interaction, 
starting from the redshift z — 3. The arrows alone the curves denote the direction of evolution. 
For both models, with the expansion of the Universe, w decreases and w' increases, and the 
w — w' diagnosis approaches to a fixed point asymptotically. For the non-interacting case, the 
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current value is (w,w') = (—0.982,-2.778 x 10~ 2 ), and the fixed point is (—1,0), which is 
same as that of ACDM model, denoted by a star symbol on the plot. Therefore, the non- 
coupling YMC model can not cross the phantom divide w = —1. For the interacting case, the 
situation is different: T = 0.31F yields a current value (w,w') = (-1.063,-5.430 x 10~ 2 ), 
and approaches to a fixed dot (—1.12,0) asymptotically. Therefore, the interaction between 
matter and DE causes a deviation from the ACDM, and give an EOF crossing the phantom 
divide w = —1. 



In the following, by using the maximum likelihood method, we will perform the best fit 
analysis on our 3-loop YMC model with the latest observational data of SNIa, CMB and LSS. 
First, we derive the constraints on the model from SNIa. Recently, the up-to-date gold sample 
of 182 SNIa data was compiled by Riess et al. [52]. It provides the apparent magnitude m(z) 
of supernovae, which is related to the luminosity distance d L (z) of supernovae through 



where M is the absolute magnitude, that can generally be considered to be the same for SNIa. 
In a flat universe the luminosity distance satisfies 



where E(z) = H(z)/H and the Hubble scale Hq 1 = 2997.9/i~ 1 Mpc. The data points of the 
latest 182 SNIa Gold dataset compiled in [52] are given in terms of the distance modulus 



4. Cosmological Constraints From SNIa, CMB, and LSS 



m(z) = M + 5\ogd L (z) +25, 



(32) 




(33) 



Hobs{zi) = m obs (zi) - M. 



(34) 



On the other hand, the theoretical distance modulus is defined as 



fith(zi) = m t h(zi) ~ M = 5 log 10 d L (zi) + 25. 



(35) 



The theoretical model parameters are determined by minimizing 



Xliv(p) = E 



i=i 



182 




(36) 
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where a is the corresponding la error, and p denotes the model parameters. In this work, we 
will determine the best fit values of corresponding model parameters (including the present 
fractional matter density f2 m0 , the Hubble constant h, and the decay rate T) in 3-loop YMC 
model. 

Next, we also consider the constraints from CMB [2] and LSS [3] observations. For the 
CMB data, we use the CMB shift parameter R, which is perhaps the most model-independent 
parameter that can be extracted from CMB data. The CMB shift parameter is given by [56] 



where the redshift of recombination z rec = 1090, which is given by WMAP5 [74]. The shift 
parameter R relates the angular diameter distance to the last scattering surface, the comoving 
size of the sound horizon at z rec and the angular scale of the first acoustic peak in CMB power 
spectrum of temperature fluctuations [56, 57]. The measured value of R has been updated 
to be R obs = 1.710 ± 0.019 from WMAP5 [74]. On the other hand, for the LSS data, we use 
the parameter A from the measurement of BAO peak in the distribution of SDSS luminous 
red galaxies, which is given by [58] 



where z b = 0.35. The SDSS BAO measurement [58] gives A obs = 0.469 (n s /0.98)-°- 35 ±0.017. 
here the scalar spectral index n s is taken to be 0.96 from the 5-year WMAP data [74]. Since 
both parameters, R and A, are independent of Hubble constant H and can be easily obtained 
from CMB and LSS observations, they provide robust constraints on DE models as useful 
complements to the SNIa data. We perform a joint analysis of latest 182 SNIa Gold dataset, 
the shift parameter R from CMB and the BAO peak measurement A from LSS to constrain 
the 3-loop YMC model. The total x 2 iS given by 




(37) 




(38) 



X — XSN + XCMB + XLSSi 



(39) 



where 



Xsn ' s given by Eq.(36), and the latter two terms are defined as 




(40) 



and 



Xlss — 



(A - A obs f 



(41) 
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The corresponding la errors are or = 0.03 and a A = 0.017, respectively. As usual, assuming 
the measurement errors are Gaussian, the likelihood function 

C oc e"* 2/2 . (42) 

The model parameters that yielding a minimal x 2 ar| d a maximal C will be favored by the 
observations. The results of our analysis are presented as the following. 

First for the non-interacting 3-loop YMC model. Utilizing the latest 182 SNIa Gold dataset 
alone, we plot the Xsn ar| d the corresponding likelihood C of 3-loop YMC in Fig. 7. It is 
found that the best-fit model parameters are fi m0 = 0.331 and h = 0.626, giving a minimal 
Xsn = 158.858. Besides, for 68% confidence level [75], the range of parameters are determined 
as Q m0 = 0.33118:8!?, flyo = 0.669l8:8io and h = 0.626 ± 0.004. In comparison, we also fit 
the ACDM model to the same SNIa data and find that the minimal Xsn a = 158.750 for the 
best-fit parameter Q m0 = 0.344 and h = 0.626. So fitting to the 182 SNIa Gold dataset alone, 
the ACDM model is slightly better than the 3-loop YMC model. 

Moreover, utilizing the combination of SNIa, CMB and LSS data, we plot the x 2 ar| d 
the corresponding likelihood C of 3-loop YMC in Fig. 8. It is found that the best-fit model 
parameters are f2 m0 = 0.289 and h = 0.634, giving a minimal Xmin = 160.317. Besides, 
for 68% confidence level [75], the range of parameters are determined as f2 m0 = 0.289t[J:oi4' 
Q y0 = 0.71llg:81| and h = 0.634 ± 0.004. As a comparison, we also fit the ACDM model 
to the same combined data. It is found that the minimal x m in,A = 162.460 for the best-fit 
parameter fi^ = 0.283 and h = 0.638. Therefore, on fitting to the combination of SNIa, 
CMB and LSS data, the non-interacting 3-loop YMC model is better than the ACDM. This 
makes the 3-loop YMC model more attractive. 

Finally, we turn to the coupled YMC model and constrain the decay rate T as the last 
model parameter. In Table 1, we list the minimal XsNmin anc ' Xmin f° r 3-loop YMC model 
with various T. For the latest SNIa data alone, the non-interaction YMC model is only slightly 
better than the coupled YMC. However, for the combination of SNIa, CMB and LSS data, 
the non-interaction YMC model is much more favored than the coupled YMC. This fact can 
be explained as follows: It is well known that most current observations, including SNIa, CMB 
and LSS, all favor the ACDM model. Therefore, to fit these observations well, a DE model 
should not deviate too far away from the ACDM model. Since the introduction of interaction 
between matter and DE will cause a deviation from the ACDM model, and a larger interaction 
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Table 1: The minimal XsN,min an d Xmin f° r 3-loop YMC model with various T. 



r 





0.12# 


0.18# 


XsN,min 


158.858 


158.877 


158.950 


Xmin 


160.317 


174.498 


192.018 



T yields a larger deviation, the coupled YMC model with a large T will not be favored by 
the observations. That is to say, to consistent with the current observations, the DE-matter 
interaction (if it exists) should be very small. 

5. Summary 

In this work, we extend our previously proposed YMC DE model to the 3-loop quantum 
corrections. This model can naturally solve the coincidence problem, and it can give, in the 
interaction form, an EOF crossing the phantom divide w = —1. Next, we study the Hubble 
parameter H, the deceleration parameter q, the statefinder diagnosis (r, s), and the w — w' 
diagnosis of 3-loop YMC model for both cases without and with interaction, and compare 
our results with other DE models. It is found that the 3-loop YMC model is much closer to 
the ACDM model than other dynamics DE models; and the introduction of the matter-DE 
interaction will make the YMC model deviating from the ACDM model. Finally, by using the 
observational data of SNIa, the shift parameter from CMB, and the BAO peak from LSS, we 
give the cosmological constraints on 3-loop YMC model. Utilizing the latest SNIa data alone, 
the best-fit model parameters of 3-loop YMC are fi m0 = 0.33llo;o2?, = 0-QQ9-omo and 
h = 0.626±0.004 (with la uncertainty); and combining SNIa, CMB and LSS data, the best-fit 
model parameters are Q m0 = 0.289±g;g^, Q y0 = 0.711±g;^| and h = 0.634 ± 0.004 (with la 
uncertainty). To fit the latest SNIa data alone, the ACDM model is slightly better than the 
3-loop YMC model; but in fitting of the combination of SNIa, CMB and LSS data, the 3-loop 
YMC model performs better than the ACDM model. This makes the 3-loop YMC model 
more attractive. In addition, the maximum likelihood analysis also shows that the interaction 
between matter and DE should be small. 

There are also other observations that would be helpful to constraint the DE models, such 
as the Chandra X-Ray observation [76], the lookback time data [77], the Gamma-Ray Bursts 
[78] and so on. In addition, it is also very interesting to constrain the YMC model by using the 
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global fitting to the full CMB and LSS data via Markov Chain Monte Carlo analysis. These 
issues deserve further investigations in the future. 
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Figure 1: Dynamical evolution of p y , p m , and p r in 3-loop YMC model without and with 
interaction. 




Figure 2: Evolution of EOS w in 3-loop YMC model without and with interaction. 
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Figure 3: Comparison of the observed H(z) in square dots with the predictions by ACDM 
model [23] and by coupled 3-loop YMC model. 




Figure 4: Evolution of deceleration parameter q(z) in 3-loop YMC model without and 
with interaction. 
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Figure 5: Evolution trajectories of stateflnder in r — s plane for 3- loop YMC model without 
and with interaction. 
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Figure 6: Evolution trajectories of w-w for 3-loop YMC model without and with inter- 
action. 
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Figure 7: The Xsn an d the corresponding likelihood £ of 3-loop YMC, where the best-fit 
parameter h = 0.626 is adopted. These results are obtained from the latest 182 SNIa 
Gold dataset. 
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Figure 8: The x 2 an d the corresponding likelihood £ of 3-loop YMC, where the best-fit 
parameter h = 0.634 is adopted. These results are obtained from the combined SNIa, 
CMB and LSS data. 
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